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Abstract 

The calculation presented in "Rotation Velocities and Radial Electric Field in the Plasma Edge" 
by W. M. Stacey [Contrib. Plasma Phys. 46 (2006)] contains an inconsistent treatment of the 
electrostatic potential. Comparing the expressions for the potential associated with the radial 
electrostatic field with that associated with the poloidal electrostatic field reveals the inconsistency. 
A field-theoretic perspective implies that the electrostatic field must vanish in a model based upon 
the physics of a neutral, conducting fluid. 
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The calculation presented in "Rotation Velocities and Radial Electric Field in the Plasma 

n 

Edge" by W. M. Stacey [1[ contains an inconsistent treatment of the electrostatic potential. 
Comparing the expressions for the potential associated with the radial electrostatic field 
with that associated with the poloidal electrostatic field reveals the inconsistency. A field- 
theoretic perspective implies that the electrostatic field must vanish in a model based upon 
the physics of a neutral, conducting fluid. 

The Georgia Tech Fusion Research Center (GT-FRC) model of Reference [lj takes the 
radial electrostatic field from the equation of motion for arbitrary species s as 

1 dp 

.... . —'ir.txr:,. „ miu - 

evaluation of that equation from experimental measurements neglects any poloidal depen- 
dence; however, when the intrinsic variation of quantities induced by the geometry is con- 
sidered JsJ, given in the large aspect ratio e = t/Rq <C 1, concentric circular R r = R flux 
surface approximation B = {0,Bg,B^} by 

B = B°/(l+£cos#) , Vg s = Vg° s /(l+ecos9) , V+. = V$ a (l + ecos9) , n s = n° s , (2) 

where A = (A) = § d9(l + e cos 9)A/2n is the flux-surface average of A. Using these values 
(r dependence implied), we find 

E r {0) = P' s °/n° s e s + V$ a B° e - ^/(l + e cos Of , (3) 

for species pressure p' s = dn s T s /dr and kT e = T e , thus E r (6) ^ E®. Expanding the denom- 
inator reveals a power series in ecosO, and as integration with respect to r does not affect 
the 9 dependency, the potential associated with the radial electrostatic field may be written 
as a cosine series 

<& Er = $° Br + ^ Er cos 9 + $| r cos 29 + . . . , (4) 

where we note the explicit absence of any sin 9 dependence. 

Turning now to the electrostatic potential appearing in the poloidal equation of motion, 
which since 1992 6| has been expressed for Eg = d&/rd9 as 

*e b = (r) [1 + $ c Eg (r) cos 9 + $% g (r) sin 9] , (5) 

where §° Eg (r) = — f a drE® ^ for a last closed flux surface at r = a, we note 
that the solution of the Fourier flux surface moments, defined by the expressions 



(A)u,c,s = § dO{l, cos6, sin 6}{1 + e cosQ)A/2ir 1 of the electron poloidal equation of mo- 
tion, T e dn e /rd6 = —en e Eg when other terms are assumed negligible at equilibrium, used 
to obtain the GT-FRC expressions relating the coefficients appearing in Eg to the extrinsic 
variations of the electron density = n^ s T e /e§° Eg , itself becomes inconsistent beyond 
zeroth order in e [7| yet will be retained herein to maintain consistency with the article com- 
mented upon. The neglect of the 0(e) terms in the system of equations indicates a solution 
valid only on the magnetic axis where e = for Rq ^ oo, yet the GT-FRC model commonly 
addresses the physics near the edge of the confinement region l|, |4| . With vanishing extrinsic 
poloidal dependence to the electron density, as in the paragraph above, the potential retains 
no explicit poloidal dependence. A sin 9 dependence to the electron density introduces a 
similar dependence in $g e but not one in $^ r , when the effect of radial dependence to the 
poloidal coefficients is neglected, as in the GT-FRC model. Thus, we conclude that the 
electrostatic potentials associated with the radial and poloidal components of the electro- 
static field evaluated from the equation of motion by the GT-FRC model are inconsistent, 

The root of the problem lies in the treatment of the radial and poloidal components of 
the electrostatic field as separate entities, when there is only one electrostatic potential from 
which both components may be calculated [8fl, stemming from the neoclassical prescription to 



relegate the defining relation for t 



moment 



0, Q, fl Q Q, Q, H, Q, 0, Q 



ie electrostatic field, Gauss' law, to a position of subsid iary 
19], a procedure fraught with difficulties 0, Q|- 



For a model of confined plasma based on the physics of a neutral, conducting fluid, where 



neutrality is understood to hold down past the scale of the differentia 



defining the continuum quantities, the electrostatic field must vanish 21[. To support an 



volume element 



electrostatic field within the medium, a neutral plasma must behave as a polarized dielectric 
along the component directions considered; however, the presence of currents in the poloidal 
plane driven by an applied toroidal electric field Q counters such hypothesis. In summary, 
a consistent treatment of the electric field within a tokamak relies on the proper application 
of all of Maxwell's laws. 
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